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Displacement of X-Ray Sources from Star Clusters in 
Starburst Galaxies 
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ABSTRACT 

We examine the spatial offsets between X-ray point sources and star clusters in 
three starburst galaxies. We find that the X-ray sources are preferentially located 
near the star clusters. Because the star clusters are very good tracers of the star 
formation activity in the galaxies, this indicates that the X-ray sources are young 
objects associated with current star formation. We find significant displacements of 
the X-ray sources from the clusters. These displacements are likely due to motion of 
the X-ray sources and indicates that they are X-ray binaries. We find that brighter 
X-ray sources preferentially occur closer to clusters. The absence of very bright sources 
at large displacements from clusters may help constrain models of the sources. 

Key words: black hole physics - galaxies: individual: M82, NGC 1569, NGC 5253 - 
galaxies: starburst - galaxies: stellar content - X-rays: galaxies 



1 INTRODUCTION 

Starburst galaxies contain prominent populations of X- 
ray sources (Fabbiano 1989 ) including unusually lumi- 
nous ones UKaaret et al. 2001llZe"zas fc Fabbiano 200211 . The 
young stellar ages within starbursts suggest that the as- 
sociated X-ray sources are likely to be young objects 
(Kilga rdeTal. 2002| . Strong X -ray variability in many 
sources further suggests those sources are accreting com- 
pact objects (Makishima ct al. 2000 Kaaret et al. 200ll 
IFabbiano et al. 2003L 

In starburst galaxies, a substantial fraction of young 
stars are found in compact, luminous star clusters 
tMeurer et al. 199511 . Indeed, Tremonti et al. (2001) and 
Harris et al. (2001) find evidence that most high mass 
stars in starbursts could form in dense clusters which 
quickly dissolve to feed the surrounding field. Consis- 
tent with this, the young compact cluster R136a con- 
tains a rich population of massive stars, evidently follow- 
ing a standard Salpeter-like upper initial mass function 
(Mass ey fc Hunter 1998| l. However, there also are indica- 
tions that some compact young massive clusters have ei- 
ther a flatter than normal upper mass function or a cut- 
off at low mass l|Sternberg 1998] |Smith fc Gallagher 2001] 
|McCrady, Gilbert, fc Gra ham 2003J . The current data are 
insufficient to distinguish between these two possibilities 
l |Mengel et al. 20 02 1 . In either case, the clusters would con- 



tain relatively more massive stars, and thus potentially more 
massive X-ray binaries. 

Star clusters in starburst regions usually are cen- 
trally concentrated, with as much as lO 6 Af0 of stars 
within a half light radius of a few parsecs (e.g. Ho 
& Filippenko 1996a, b). Stellar encounters in such dense 
clusters may enhance the production of X-ray binaries 
l |Portegie s Zwart et al. 1 999^. Interact ions of binaries in 
clusters ( |PlimnByT^^igurdsson 1991| and binary recoil 
following supernovae leading to compact object forma- 
tion jNclcmans, Tauris, fc van den Heuvel 1999T > can lead to 
ejection of X-ray binaries from their point of origin. There- 
fore, if the X-ray sources are X-ray binaries, then they may 
be expected to be spread over a larger spatial area than the 
star clusters. 

Here, we study the spatial offsets between X-ray point 
sources and star clusters in three starburst galaxies with 
dense clusters and deep, high-resolution X-ray observations. 
We describe our sample, the observations, and our analysis 
in § 2, the results in § 3, and discuss the implications in § 4. 



2 DATA 

Our selected galaxies are: M82, NGC 1569, and NGC 5253. 
All three are relatively nearby, within 4 Mpc. Hence, rela- 
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Figure 1. False color infrared image of M82. The NIC2 data were used for the inner parts of the mosaic and the NIC3 for the remainder. 
Candidate star clusters are marked as diamonds. X-ray sources are marked as X's. There are 8 sources in the sample which do not appear 
on the figure. North is up, and the axes are R.A. and DEC in J2000 coordinates. 



tively dim X-ray sources can be detected with modest expo- 
sures. NGC 5253 is a dwarf galaxy which has been undergo- 
ing a starburst for no more than 10 Myr l|Beck et al. 1996")l 
and the young star clusters have ages between 1 Myr and 
8 Myr (ICalzetti et al. 19971 ITremonti et al. 200111 . M82 is 
the prototypical starburst galaxy. The nuclear star burst 
has been active for roughly 10 Myr (Satyapal et al. 1997| l 
and may consist of two events with ages of 5 Myr and 
10 Myr ( |McCrady, Gilbert, fc Graham 20031 . There is also 
a more spatially extended "fossil" starburst which occurred 
0.6-1 Gyr ago (|de~Grijs, O'Connell, fc Gallagher 2001| 
|de Grijs, Bastian, fc Lamers 2003) . NGC 1569 is a dwarf ir- 
regular galaxy which has been in a starburst phase for 10- 
20 Myr (IHunter et al. 200011 . 

2.1 Star clusters 

For each galaxy, we require the positions of the star clusters. 
For NGC 1569, we used the list of star clusters compiled by 
Hunter et al. (2000) based on HST observations in the opti- 
cal band. Because this galaxy has a relatively low gas/dust 
content and therefore low internal extinction, the optical 
data should provide a complete cluster list. 

M82 is a larger galaxy and has high obscuration 
near its core, where most of the star clusters are located 
(McLcod ct al. 1993). Comparison of optical and IR images 
shows that many clusters are missed in the optical due to 
the high obscuration. We chose to derive a cluster list from 
IR observations obtained with NICMOS on HST. The ob- 
servations and reduction of data to images are described in 
Alonso-Herrero et al. (2003). We used a sliding cell algorithm 
to detect compact sources within the images from the NIC2 
camera using the F160W filter with a central wavelength 
1.6 fim and covering a wavelength range of 1.4-1.8 /^m, and 
the NIC3 camera using the F166N filter with a central wave- 
length 1.66 (jm and a 1% bandpass (images in the F160W 
filter were not available with NIC3). The NIC2 images are 



higher resolution (75 mas pixels) and cover the inner regions 
of the galaxy. The NIC3 images extend the imaged field at 
lower resolution (200 mas pixels). 

We found clusters in each image. Clusters found in mul- 
tiple images were used to align the relative positions of the 
various NICMOS images. We determined the absolute as- 
trometry of the NICMOS images using the positions of four 
clusters detected in the 2MASS survey which lie well out- 
side the crowded central region of the galaxy. After applying 
one global shift to the NICMOS coordinates, the positions 
of all four clusters agree within 0.5" with the 2MASS source 
positions. A mosaic of the NICMOS data is shown in Fig.0 

For NGC 5253, we also used NICMOS observations 
taken using the NIC2 camera and the F160W filter. The 
details of these observations, the analysis, and the cluster 
list will be presented elsewhere (Alonso-Herrero et al., in 
preparation). The NICMOS field of view is 19" x 19" and 
covers the central region of the galaxy. The astrometry was 
set using the accurate position for the dominant 1.3 cm radio 
source l |Turner, Beck, fc Ho 200*01 which is coincident with 
one of the IR clusters. Due to the small field of view, the 
roll angle uncertainty from the HST aspect solution leads to 
only small uncertainties in the cluster positions. 

2.2 X-ray sources 

We extracted the longest available observation of each 
galaxy from the Chandra X-Ray Observatory data archive. 
All observations employed the Advanced CCD Imaging 
Spectrometer and the High-Resolution Mirror Assembly. 
The data for M82 JcTriffiths et al."2000l and NGC 1569 
l |Martin, Kobulnicky, Heckman 2002^ have been published 
previously. We reanalyzed the data to have consistent re- 
sults for all three galaxies. 

For each galaxy, we extracted an image of X-rays in 
the 0.3-8 keV band. We calculated an exposure map for 
an assumed powerlaw spectrum with interstellar absorption. 



© 2004 RAS, MNRAS 000,[I]|!>] 



X-Ray Sources and Star Clusters 3 



Table 1. Statistics of the Displacement Distributions. 



Galaxy 


X-ray 
Sources 


Clusters 


Displacement (pc) 
Average Median 


M82 

NGC 1569 
NGC 5253 


42 
11 
10 


50 
58 
13 


201 ± 37 
218 ± 95 
239 ± 106 


118±ii 

34±r 

90+ 2 * 

— o2 



The Table contains for each galaxy: the galaxy name, the number 
of X-ray sources within 1 kpc of a star cluster, the number of star 
clusters, and the average and median of the displacements of X- 
ray sources from star clusters. 



The powerlaw photon index was fixed to 1.8 in all cases, and 
the absorption column density set to the Galactic line-of- 
sight value for each galaxy. We detected sources using the 
routine wavdetect in CIAO version 2.3 requiring a source 
significance of 3<r. We compared the resulting source list with 
optical images from the POSS2 to remove foreground stars 
and with the NED extragalactic object database to remove 
background AGN. We used X-ray sources with either AGN 
or optical star counterparts (from the USNO B catalog) to 
check the astrometry. The absolute X-ray source positions 
should be accurate to better than 1". 

Because spectral fitting is not feasible except for the few 
brightest sources, we calculated an X-ray flux using the same 
spectrum assumed for the exposure map and correcting only 
for interstellar absorption within our Galaxy. If the sources 
have intrinsic absorption or there is absorption within the 
host galaxy, the true fluxes may be higher. We calculated 
a luminosity for each X-ray source assuming that it is at 
the distance to the corresponding galaxy which we take to 
be 3.6 Mpc for M82, 2.2 Mpc for NGC 1569, and 3.3 Mpc 
for NGC 5253 HGibson et al. 2 000). Note that we calculate 
luminosities assuming isotropic emission. If the X-rays are 
beamed, then the luminosities would be lower. The lumi- 
nosities are calculated for the 0.3-8 keV band. The weakest 
detected sources are 3 x 10 36 erg s -1 for M82, 5 x 10 35 erg s _1 
for NGC 1569, and 1 x 10 36 ergs" 1 for NGC 5253. 

While some of the X-ray sources may be supernova 
remnants (SNRs) or young supernovae, the fraction of such 
sources is likely to be small. Martin et al. (2002) identify only 
one Chandra X-ray source in NGC 1569 as a possible SNR. 
Griffiths et al. (2000) report three close coincidences (within 
1") of Chandra X-ray sources with radio SNRs in M82. How- 
ever, the sources have hard X-ray spectra rather the soft 
spectra which are found for SNRs fPrest wicFet al. 200311 
and young supernovae (Kaaret 2001). 

3 RESULTS 

To investigate the relation between star clusters and X-ray 
sources, we found the nearest star cluster to each X-ray 
source. For X-ray sources born in clusters, this provides a 
lower bound on the displacement of the source from its par- 
ent cluster. The uncertainty in the relative alignment of the 
X-ray source and star cluster positions is 1" or less, cor- 
responding to uncertainties in the displacements of about 
10 pc. Displacements smaller than 10 pc have been set to 
10 pc. Displacements smaller than about 20 pc should be 
considered as upper bounds. We examine only sources within 
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Figure 2. Cumulative distribution of X-ray source displace- 
ments for sources in the luminosity range lO^ergs" 1 < L x < 
10 38 ergs- 1 . The solid line is for M82, the dotted line for NGC 
1569, and the dashed line for NGC 5253. Each histogram is nor- 
malized by the number of X-ray sources within 1 kpc of a star 
cluster for that galaxy. 



1 kpc of a star cluster. The distributions of the spatial dis- 
placements are shown in Fig. |2] and some statistics of the 
distributions are shown in Tabled 

Our cluster lists contains all of the most luminous clus- 
ters for each galaxy, but may be incomplete at low luminosi- 
ties. To determine if our results are sensitive to the com- 
pleteness of the sample, we varied the detection threshold 
for the clusters in M82 and found the distributions of the 
spatial displacements of X-ray sources from the clusters for 
the various cluster lists. We found no significant changes 
in the distribution for reasonable changes in the threshold. 
This may be because the lower luminosity clusters tend to 
be located near brighter clusters already in the list. Thus, 
the displacements are decreased only slightly. 

To determine if the clustering of the X-ray sources near 
the star clusters is statistically significant, we generated ran- 
dom sets of uniformly distributed sources and found the 
spatial displacements from the star clusters following the 
same procedures used for the actual X-ray sources. The dis- 
placement distributions of the X-ray sources are inconsistent 
with that expected for a uniform distribution of sources. The 
probabilities of chance occurrence from a uniform distribu- 
tion are 2 x 10 -6 for M82, 5 x 10" 5 for NGC 1569, and 

2 x 10 -4 for NGC 5253. For each galaxy, the average dis- 
placement of X-ray sources from star clusters is significantly 
smaller than for the random source distribution, indicating 
that the X-ray sources are preferentially located near the 
star clusters. 

The distributions of the spatial displacements are sim- 
ilar in all three galaxies, see Fig. [5] We compared the vari- 
ous distributions using a Kolmogorov-Smirnov (KS) test and 
also by calculating the average source displacement. We ap- 
plied both tests to both the full data sample and also re- 
stricting the source luminosities to the range 10 36 ergs _1 < 
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Figure 3. X-ray source luminosity versus displacement from 
nearest star cluster. Sources from M82 are shown as diamonds, 
NGC 1569 as triangles, and NGC 5253 as squares. The luminosi- 
ties are for the 0.3—8 keV band. The dashed line represents the 
equation L x = (1 X 10 41 erg s~ 1 )/(d/pc). 



L x < 10 38 ergs _1 . In all cases, the distributions of the spa- 
tial displacements for the X-ray sources in the three galaxies 
are consistent with being drawn from the same distribution. 
The median displacement for NGC 1569 appears somewhat 
smaller than for the other galaxies, but the difference is not 
statistically significant. The 90% confidence level error in- 
tervals for the median overlap for all three galaxies. 

Fig.[3]shows the luminosity of each X-ray source plotted 
versus its spatial displacement from the nearest star cluster. 
The striking feature of the plot is that there are no high 
luminosity sources at large displacements from star clusters. 
Furthermore, there is an apparent trend of decreasing lu- 
minosity with increasing displacement from the nearest star 
cluster. 

To evaluate if this trend is statistically significant, we 
performed a 'boot-strap' analysis. We used the set of dis- 
placements and luminosities shown in Fig. [3] and then ran- 
domly re-arranged the pairings. This corresponds to the null 
hypothesis in which luminosity and displacement are unre- 
lated. In 10 7 trials, we found 1920 cases in which no sources 
lay in the same region of the luminosity-displacement plot 
which is empty in Fig. [3] This corresponds to a chance prob- 
ability of occurrence of 1.9 x 10 -4 . 

For M82 there is also an absence of dim sources 
at small displacements. This is likely due to the high 
level of diffuse X-ray emission in the central 200 pc of 
M82 ( Griffiths ct al. 2000 1 which precludes detection of dim 
sources in the central regions of M82. Dim sources are de- 
tected very near star clusters in NGC 1569. 



4 DISCUSSION 

We have shown that X-ray sources in these three starburst 
galaxies are preferentially located near star clusters. Because 



the star clusters are very good tracers of current star forma- 
tion activity in the galaxies, this confirms that the X-ray 
binaries are young objects associated with current star for- 
mation. 

We also found significant displacements of the X-ray 
sources from the clusters. Because X-ray binaries, unlike 
other bright X-ray sources such as young supernovae, can 
exhibit high velocities, this suggests that much of the 
X-ray source population consists of X-ray binaries. Several 
mechanisms can put a binary in motion. For neutron star 
binaries, a 'kick' due to an asymmetric explosion in the 
formation of the neutron star can lead to high velocities 
l |Lyne fc Lorimer 1994} . Even in the absence of 'kicks' from 
supernova explosions, momentum conservation following a 
symmetric ejection of matter in the formation of the neutron 
star or black hole in a binary with a high mass companion 
can produce a runaway speed of ~ 50kms _ ; the ejected 
matter continues to move with the instantaneous orbital ve- 
locity of the compact object at the moment of ejection and 
the binary must move in the opposite direction to conserve 
momentum l|Nelem ans, Tauris, fc van den Heuvel 19991 
Ivan den Heuvel"etwu f ^b00ll . For binaries in clusters, inter- 
actions with other stars and binaries in the cluster can eject 
the binary from the cluster (jPhinney fc Sigurdsson 1991| 



Kulk arni, Hut, fc McMillian 1993| " 
Slgurd^sonT^HernquisT^^93^ In young, dense star 



clusters such interactions can occur on time scales 
of a few Myr ( Portegie s Zwart et al. 19 99 1 . Objects 
ejected via dynamical interactions tend to escape 
with close to the minimum energy needed to escape 
<Joshi, Nave, fc Rasio 2001) . The runaway velocities (at 
infinity) should be of the same magnitude as the stellar 
velocity dispersions of the clusters, which are typically 
10 — 15kms -1 in these galaxies (|Smith fc Gallagher 2001| 
|McCrady, Gilbert, fc Graham 2003| . The displacements 
we observe in the starburst galaxies are likely due to 
motion of the X-ray sources caused by one or more of these 
mechanisms. 

Furthermore, we found that there is an absence of bright 
X-ray sources with large displacements. This suggests that 
there is some correlation between the maximum possible 
brightness of an X-ray source and its motion. This corre- 
lation appears to hold only for (isotropic equivalent) lumi- 
nosities above 10 38 ergs -1 . The excluded region appears to 
be bounded by a linear relation between (isotropic equiv- 
alent) source luminosity Lx and source displacement from 
the nearest star cluster d, L x < (1 x 10 41 ergs _1 )/(d/pc). 

In discussing this correlation, we first consider the case 
where the X-ray sources emit isotropically. In this case, 
the systems producing such high luminosities likely con- 
tain black holes accreting via Roche lobe overflow because 
such high luminosities would be difficult to achieve in a 
wind accretor due to the low efficiency of wind capture and 
black holes and needed to not violate the Eddington limit 
jBlondin, Stevens, fc Kallman 1991| IPetterson 19781 . 

If the X-ray sources are ejected from the star clusters 
with speeds which are roughly independent of mass, then 
the inverse correlation between maximum X-ray source lumi- 
nosity and displacement from the nearest star cluster would 
arise if the source lifetime varies inversely with luminosity. 
An upper bound on the source lifetime can be obtained from 
the time required to fully accrete the stellar companion. For 
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a companion mass M and an efficiency for the conversion of 
mass lost by the companion to luminosity of 77, the source 
lifetime must be T ^ r/Mc 2 /L where L is the average lumi- 
nosity and c is the speed of light. For sources traveling with 
a speed v perpendicular to the line of sight, the displace- 
ment from the point of origin will then be d ^ vrjMc? / L. If 
the companion mass is independent of the compact object 
mass, then this would reproduce the required dependence of 
source lifetime on luminosity. 

Given a typical runaway velocity v ~ lOkms -1 , we 
must have r\M ~ 0.2Mq to match the line plotted in Fig. [3] 
which bounds the region where X-ray sources are found. If 
Roche lobe overflow is occurring, then accretion onto the 
compact object may be efficient with little mass loss giv- 
ing 77 ~ 0.1. In this case, the companion mass would be 
M ~ 2Mq. Such intermediate mass companions could be 
captured via dynamical interactions in the cluster. However, 
the capture must be directly into a Roche-lobe filling orbit or 
the binary must hardened into a Roche-lobe filling orbit via 
successive interactions in order to begin accretion promptly, 
since the evolutionary time scale of the companion is long. 
A better understanding of the IMF and the dynamical in- 
teractions within the clusters is needed to determine if this 
scenario is viable. Even with accretion via Roche lobe over- 
flow, the efficiency r\ may be less than 0.1 since outflows 
are often observed in X-ray binaries. Super-Eddington mass 
transfer would also produce r\ < 0.1. If r\ < 0.1, then a higher 
companion mass may be compatible with the data shown in 
Fig.H 

If the X-ray sources have high-mass companions and 
maximum speeds near ~ 50 km s~ 1 , then the absence of high 
luminosity source at large displacements implies a limit on 
the X-ray emitting lifetime of the sources. We find no sources 
at luminosities above 10 38 ergs _1 at displacements larger 
than 200 pc. At SOkms^ 1 , this would imply that the life- 
times of these luminous sources must be less than 4 Myr, 
corresponding to very massive stars. An alternative is that 
the luminosity of the sources decreases with age. This would 
require an evolutionary path for binaries which produces an 
accretion rate which decreases with age. 

The X-ray binaries may also be beamed 
{King et al. 2001| |Kording et al. 2002| IKaaret et al. 200311 . 
If the X-ray sources are high-mass systems with high 
velocities, then the observed correlation would imply that 
beaming only occurs when the binaries are quite young. 
King et al. (2001) suggest that the ULXs in starburst 
galaxies are high mass X-ray binaries with beamed X-ray 
emission in a phase of thermal-timescale mass transfer. 
The delay between the formation of the black hole (and, 
presumably, the start of the binary's motion away from 
it point of origin) and the onset of the thermal-timescale 
mass transfer phase depends on the stellar evolution of 
the companion. The delay could be ~ 20 Myr for a 9Mq 
companion, which would imply that very bright X-ray 
sources should be visible out to 1 kpc. The data appear 
inconsistent with this, unless highly beamed X-ray emission 
occurs in the thermal-timescale mass transfer phase only 
for very massive companions. Large displacements of high 
flux sources, inconsistent with the data, also appear allowed 
in the relativistic beamed model of Kording et al. (2002). 
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